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Abstract 

We present a non-linear elastic model of a coherent transition with discontinuous 
volume change in an isotropic solid. The model reproduces the anomalous thermody- 
namics typical of coherent equilibrium including intrinsic hysteresis (for a pressure 
driven experiment) and a negative bulk modulus. The novelty of the model is that 
the statistical mechanics solution can be easily worked out. We find that coherency 
leads to an infinite-range density-density interaction, which drives classical critical 
behavior. The pressure width of the hysteresis loop shrinks with increasing tem- 
perature, ending at a critical point at a temperature related to the shear modulus. 
The bulk modulus softens with a 1/2 exponent at the transition even far from the 
critical point. Many well known features of the phase diagram of Ce and related 
systems are explained by the model. 
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1 Introduction 



Room temperature 7 cerium shows a well known transition to isostructural a 
cerium at ~ 0.8GPa with a 17% reduction of the volume. Similar (though not 
necessarily isostructural) transitions occur in other /-electron systems usually 
associated with mixed- valence behavior [1,2]. Volume-collapse transitions were 
also predicted to occur in colloidal systems [3] and close to first order electronic 
transitions when mesoscopic inhomogeneous electronic states are frustrated by 
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the long-range Coulomb interaction[4,5]. In Ce and SmS the volume collapse 
transition ends at a critical point, making of these and related systems the 
solid state analog of a van der Waals liquid-gas system. 

The thermodynamics of coherent structural transitions have been shown to 
be highly anomalous[6,7,8,9,10,ll,12,13,14,15,16]. In the mixed phase the sys- 
tem behaves non extensively and violates usual stability criteria for extensive 
systems [7, 14, 15]. In particular, it possesses negative compressibility [14, 15]. For 
an experiment in which the pressure acts as the control variable usual nucle- 
ation is forbidden within a finite pressure window. This pressure window im- 
plies an intrinsic hysteresis loop as pointed out by Roytburd[8,9,14,15,17] and 
discussed also in metal-hydrogen systems [9, 12]. Outside the pressure window 
the transition can occur, but only irreversibly with an intrinsic jump of the 
thermodynamic potential. 

In this work we present a simple non-linear elastic model describing a coherent 
volume-collapse transition in an isotropic solid. The assumed elastic coherence 
leads to a distance-independent density-density interaction analogous to com- 
positional effects in alloys[6, 7,8,9, 10] and hydrogen density-density interaction 
in metals[ll]. Mean-field (MF) theory becomes exact and the statistical me- 
chanics solution of the model can be easily worked out, which constitutes 
the main novelty of the model. The solution reproduces the above mentioned 
anomalous thermodynamics. In addition some novel results are found. The ex- 
tent of the pressure hysteresis reduces as a function of temperature, vanishing 
at a critical point. The temperature of the critical point is directly related to 
the shear modulus of the material. We show that critical exponents are MF 
like, and this remains true even in the presence of short-range interactions, 
providing a straightforward explanation to old experiments in mixed-valence 
systems [2]. We discuss under what conditions the transition occurs with a 
softening of the bulk modulus and discuss its critical behavior. 



2 Model 



In the present model, the order parameter (OP) is taken to be the dilation 
strain e = J2i=i e a where are the components of the infinitesimal strain 
tensor [18] e and d is the dimensionality of the system. The dilation strain 
describes changes in the local specific volume with respect to a reference state, 
and it is defined in such a way that if Vq is the volume of an original piece of 
material (at e = 0), then the volume for e ^ is V — V (l + e). 

Bistability in the volume is described by a coarse-grained energy functional 
Hi. If we assume that upon collapse, the system passes from a specific volume 
Vq to some lower value Vq , then an appropriate description of the transition 
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will contain a local term of the form 



H 




ae 



) 



(1) 



Here the reference state (for which e = 0) was taken to correspond to a specific 
volume v = (vq + Vq)/2, and the relative volume change at the transition (at 
zero temperature) is given by 



The linear term in (1) is included in order to have the transition occurring 
at a finite pressure. Hi alone describes an incoherent volume collapse tran- 
sition with nominal volume change Av/v and at a nominal pressure P = a. 
Coherence effects will change this simple picture. 

In a coherent transition non-uniform dilation strains are accompanied by shear 
strains. The latter give an additional elastic contribution H h to the total en- 
ergy. To lowest order we take this contribution to be harmonic, and for the 
assumed isotropic case we have: 



where \i is the shear modulus and the traceless strain tensor e^- = — b^efd 
has been defined. 

The coefficients in (1) and (3) are assumed to be temperature and pressure 
independent, i.e., Eq. (1) has to be considered as an energy functional (rather 
than a Ginzburg-Landau free energy). The phase transition will be derived 
from the present energy functional by its introduction into a partition function 
and a full statistical mechanics calculation. Note also that we are assuming 
that the two phases have the same bulk and shear modulus. 



3 Zero temperature theory 

In a coherent structural phase transition the strain tensor must fulfill the St. 
Venant compatibility constraints Vx (Vxe) 1 " = [8,9,13,14,15,17,19,20,21,22]. 
This leads to constrains between the shear strains and the OP. We take peri- 
odic boundary conditions and work with the Fourier transform of the strain 
tensor, e(k). 

The k = and k ^ components have to be treated separately [23]. Indeed a 
uniform strain satisfies the St. Venant compatibility constraints trivially and 






(3) 
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Fig. 1. (a) ^-dependent common tangent construction. The cases 4> = —0.5 and 0.2 
are illustrated (m = 1). The triangles indicate the minimum value of h attainable 
for the corresponding values of <f>. (b) Analytical results for h((f)) as a function of 
4> for 7 = and different values of m, as indicated. For m > 2 the solution is 

h(4>) = h (4>). 

therefore it is unconstrained. For k ^ the shear strains can be eliminated 
in favor of e by minimizing Hh with respect to the shear strains as shown 
in the Appendix. The shear strains are said to "accommodate" to the order 
parameter so they will be also referred to as the "accommodation strains" . As 
shown in the Appendix this allows us to express Eq. (3) as a function of the 
OP alone: 

^=(l-i)/i^Ee(k)e(-k). (4) 

Here L is the linear dimension of the system and e(k) are the Fourier compo- 
nents of e. A crucial point is that, since the k = term is unconstrained, it is 
excluded from the sum. We can formally extend the sum to all k by explicitly 
subtracting the missing term. One obtains in real space: 

H h =(l--^j d d x[e(x) - e] 2 (5) 

with e = / d d xe/L d . A similar result has been obtained in d = 2 in the context 
of the dislocation mediated melting problem[23]. 

Eq. (5) manifestly satisfies the Bitter-Crum theorem of coherent mixtures 
which states that the elastic energy is independent of the shape and location 
of the phases[10]. 

We define now = e^Jc/b, m = 4(d — l)/i/(db) and 7 = 2n/b. Up to an 
additive constant the coarse grained energy density h [from Eqs. (1) and (5)] 
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in units of the height of the double well barrier [b 2 /(Ac)] reads: 



h(<f>, 0) = ho{</>) + 7l V0| 2 + m(0 - 0) 2 (6) 

with ho((f)) = (1 — <f) 2 ) 2 being the double well in the dimensionless volume 0. 
The linear term in the OP has been eliminated by redefining the origin of 
pressure, i.e., the transition pressure is now zero. The parameter 7 fixes the 
width of the interface between the low and high volume phases in a mixed 
state. The energy density Eq. (6) couples the values of at different spatial 
positions only through 0, namely it corresponds effectively to a MF or an 
infinite range interaction. We stress however, that this is the exact energy 
of the model, the mean-field contribution coming from the harmonic elastic 
energy of the accommodation strains, Eq. (5). 

We first discuss the 7 = and zero temperature case. Mixed phase solutions 
for a constrained total volume obey the familiar common tangent construction: 
dh((f),<f>)/ d<f)\ 4>=( f ) _ = dh(<f>, 0)/<90l0=</, + = [/i(0_, 0)-/i(0+, 0)]/(0_-0+) where 
<p± = ±y (1 — m/2) and the + (— ) corresponds to the OP of the expanded 
(collapsed) phase. The only difference with a conventional phase separation 
computation is that the construction is 0-dependent because h depends on 
the global variable [see fig. 1(a)]. This dependence makes the total energy of 
the mixture non additive. In fig. 1(b) we report the resulting average energy 
density h(4>). 

For m > 2 the mixed solution does not exist and the system remains always 
uniform. For m < 4/3 the mixed state solutions produce a shift of the spinodal 
points of the system [fig. 1(b)] whereas for m > 4/3 the spinodal points 
coincide with the spinodal points of ho. 

The dimensionless pressure p can be calculated as p = —dh((f))/d(f). The phys- 
ical pressure P is obtained from p as P = p^b 3 / (16c) + a. The nominal 
transition pressure is p = (i.e., P = a). In the mixed phase the free en- 
ergy has negative curvature so the compressibility is negative [fig. 1(b)]. This 
is also true for uniform phases for m > 4/3 with volume between the spin- 
odal points. This thermodynamic state is stable at constant volume with the 
pinned boundary conditions of Ref. [24,25]. We mention that this violation of 
usual stability criteria is common in other contexts where long-range forces 
appear [2 6]. 



4 Finite temperature theory 

For the model with 7 = 0, it is straightforward to obtain the exact partition 
function and the exact — p dependence at finite temperatures. Let us rewrite 
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for convenience the energy density h in the form 



h{<j>,<j>) = h {4>) +m(0-0) 2 + P 0, (7) 

where we set 7 = and introduce explicitly the pressure term. 

To solve the model we have to consider the minimum volume that can fluctuate 
independently in the system. This volume will be noted £ d , with £ an atomic 
length scale. The energy is written in discretized form as a sum over N cells 
of volume £ d with L d = N£ d : 

h 2 N 

The partition function can be written as Z = z N with 

z{4>) = J d4>e- h ^'\ 

where t is a dimensionless temperature defined as: t = ksTAc/ \b 2 ^ d ). The 
value of must be computed self-consistently as: 

4> = -L / d^e-^^'K (8) 

z{<j>) J 



In fig. 2(a) we show the isotherms in the 0-p plane for m = 6 for different values 
of the dimensionless temperature t. Controlling the total volume the system 
would follow the dashed lines. For t = we also show the hysteresis loop 
(full line) that the system would follow in a pressure controlled experiment. 
Starting from the high volume phase at low pressure and increasing p the 
system remains in this phase until the spinodal pressure , i.e. until the limit 
of stability of the phase where the volume suddenly collapses. Reducing the 
pressure the system suddenly expands at pK A conventional (incoherent) first- 
order phase transition driven by pressure at a very low rate (so that nucleation 
occurs) would follow a vertical line at p — 0. This is not possible here because 
the infinite range interaction makes the critical drop radius[18] diverge, and 
the nucleation energy-barrier scale as the volume of the system, a well known 
fact in structural transformations [16]. 

The pressure width (Ap) and volume width ( A0) of the hysteresis loop as com- 
puted from Eq. (8) are reduced at finite temperatures up to the point where 
they vanish, which determines the critical point coordinates (p c = 0,t c (m)). 
The pressure width as a function of temperature determines a wedge of forbid- 
den nucleation [shown in fig. 2(b) for different values of m] around the critical 
pressure p c = 0. 
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Because of the mean-field character of the model the isotherms are analytic 
functions when seen as pressure vs. volume plots [rather than volume vs. 
pressure as in Fig. 2 (a)]. Therefore close to the spinodal, inverting the analytic 
behavior, one finds that the volume behaves as [p— pi- 1 -] 1 / 2 . The compressibility 
(oc d(f)/dp) diverges at the transition even far from the critical point with a 
1/2 exponent. 

The solution in the case of finite 7 and m is given implicitly by = (p (p + 
m4>) where 0oG°) is the equation of state of the same model with m — 0. 
This corresponds to the well studied 4 -model which is in the universality 
class of the Ising model[18]. We can provide here only a description of the 
results that are obtained, a full detailed solution will be presented elsewhere. 
If m > 4/3 the qualitative description of the forbidden nucleation region are 
completely equivalent to the 7 = case, i.e., the previously discussed pressure- 
temperature phase diagram and in particular the description of the forbidden 
nucleation wedge are not qualitatively affected. For m < 4/3 and sufficiently 
low temperatures, by driving the system through the transition, the limit 
of forbidden nucleation is reached before the stability limit. Therefore, there 
exists a pressure region, outside the window of forbidden nucleation, where 
the system can still remain in a metastable state until reaching the spinodal 
pressure, and a standard nucleation and growth process may occur due to the 
surface energy contribution of the gradient term. However, it is possible to 
show that this nucleation and growth region never reaches the critical point. 
From this we can conclude that classical critical exponents are obtained also 




P P 



Fig. 2. (a) (f> vs. p for m = 6 and different values of the dimensionless tempera- 
ture t, as indicated, (b) t-p phase diagram for 7 = and different m values, as 
indicated. The interior of the wedge (for each value of m) is a region in which nu- 
cleation is forbidden. The crosses are the experimental width of the hysteresis loop 
for Smo.87Gdo.13S [27] rendered dimensionless with theoretical estimates of b and c. 
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for m < 4/3. Thus we find that critical exponents take the mean- field values 
for both 7 = and 7 > . This is also in agreement with the d = 2 theoretical 
analysis of Ref. [23] . 



5 Comparison with experiments 

A wedge where nucleation does not occur is observed in Ce[l,2], SmS[27] and 
related systems. Interestingly it is also observed in a volume-collapse transition 
in amorphous ice [28]. However without detailed studies it is difficult to be sure 
that it does not arise from the usual hysteresis due to a finite driving rate, as 
occurs in other first order transitions. Other indicators, however, point more 
strongly to the applicability of our theory. 

In agreement with the present theory MF critical exponents have been ob- 
served in Ce systems [2]. The connection with long-range strain effects has 
long been suspected[2] but, to the best of our knowledge, has not received a 
microscopic explanation in mixed-valence systems. 

In contrast with usual first-order phase transitions, where hysteresis is not 
intrinsic, one can define here a critical exponent for the pressure width of the 
hysteresis loop: Ap oc \t — t c \ 3 ^ 2 . In particular the curvature of the boundaries 
close to the critical point must be convex. This is in agreement with the 
experimental data on Gd doped SmS[27] shown in Fig 2(b) supporting our 
conclusion that Ap is an intrinsic property of the system. 

Also in contrast with usual first-order transitions, critical behavior (although 
classical like) can be easily observed far from the critical point. The diver- 
gence of the compressibility discussed above implies that the bulk modulus 
vanishes as \p — p^\ 1 / 2 and jumps discontinuously to a finite value for p > p\ 
Remarkably, after this work was completed and a preliminary version of it was 
posted[29] this prediction was experimentally verified in Ce[30]. 

Decreasing the pressure, the jump and the square root singularity change sides 
and occur at pK Similar behavior is found as a function of temperature. Also 
an anomaly measured in YbInCu4 on cooling resembles our prediction when 
plotted in an appropriate scale (see inset of fig. 4 in Ref. [31]). This also 
explains the appearance of precursor effects which have been very puzzling 
in the mixed-valence literature and have led some authors to interpret the 
transition as second order [31]. 

A divergent compressibility implies the unusual property that the Poisson's 
ratio will become negative and approach — 1. A negative Poisson's ratio is 
indeed observed in the SmS system and is generally associated with mixed- 
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valence behavior [32]. Since the unusual elastic properties we find (including 
negative compressibility) occur in extensive temperature and volume ranges 
which can additionally be tuned by alloying and external fields we conclude 
that mixed-valence systems are ideal candidates for applications like induced 
hardening and giant damping (c.f. Ref. [25]). 

In mixed- valence systems the double well is attributed to the anomalous con- 
tribution of /-electrons to the crystal binding [33, 34]. For a realistic description 
it may be necessary to take the parameters in Eq. (1) to depend on temper- 
ature to take into account temperature effects of the /-electrons. In fact a 
trivial linear and quadratic temperature dependence of a fits our symmetric 
phase diagram to the asymmetric forms of Ce[l,2] and SmS[2,27]. 

The temperature dependence of a does not affect the critical temperature, 
but b does. Indeed in current explanations of the critical point in /-electron 
systems this electronic effect is crucial whereas the shear rigidity does not play 
a relevant role [33, 34]. 

To check the relevance of the temperature scale set by the shear modulus 
in determining real values of T c we take the opposite point of view, i.e., we 
neglect the temperature dependence of the coefficients in Eq. (1) and take into 
account the effect of the shear rigidity alone. 

We can eliminate the unknown coefficients a and b in favor of the relative vol- 
ume Av(T)/v and pressure jump AP(T) at a given temperature T. Av(T)/v 
is defined at the pseudo equilibrium pressure (practically the average between 
P^ and P- 1 ). Restoring dimensions [ksT = t^ d b 2 /(4c)} the critical temperature 
for m ^> 4/3 is given by: 



2 

+ T5T> (9) 

k B 5T' = CAP(T)^f±£ d , k B 5T ee (^v 1 ) 2 and C ~ 0.137. 

This relation between \x and T c is rooted in the fact that the critical point along 
the critical isochore can be seen as an order-disorder transition governed by 
the rigidity of the lattice parametrized by //, in strong contrast with current 
theories of the critical point in Ce[33,34]. A similar relation between \x and T c 
has been obtained for metal hydridesfll]. 

For the 7 = case, £ d should be assimilated to the atomic volume. In the 
case 7 > 0, this volume can be greater depending on the strength of the short 
distance interactions parametrized by 7. In any case we expect it to be of 
the order of a few formula unit volumes Vf. To make a rough estimate we 
took £ d = 2vf and we obtain good agreement with the experimental values 



m m 5T + 8T 



'ST + 8T> 
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Table 1 

Parameters at T = 293K used to estimate T c [K]. AP and \i are given in GPa. fi 
is an average value. Experimental data for Smi-^Gd^S after Ref. [27] except for \i 
which is taken at x = 0[35]. We caution that the experimental values of T c are not 
well established and vary as much as 10% for different authors. 





A* 


Av/v 


AP 


v f [A 3 ] 


rpExp. 


rpTheory 


m 


Ce[l,2] 


13 


0.17 


0.2 


31.3 


550 


530 


5 


SmS 


40 


0.105 


0.2 


53.3 


700 


702 


12 


Smo.87Gdo.13S 




0.075 


0.2 


50.5 


543 


508 
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(Table 1). 

Since the relation between £ d and Vf goes beyond our model we cannot exclude 
temperature effects on the coefficients of Eq. (1) contributing to T c [33,34]; 
however, the fact that we obtain the right order of magnitude with a £ d of 
the order of Vf implies that the scale set by the shear modulus can not be 
neglected in T c . 



6 Conclusions 

In conclusion we have presented a simple solvable model of an isostructural 
volume collapse transition, the solid-state analog of the van der Waals tran- 
sition. The model reproduces the anomalous thermodynamics found in other 
approaches[6,7,8,9,10,ll,12,13,14,15,16]. In particular, the transition is intrin- 
sically irreversible with a discontinuity of the thermodynamic potential. The 
main novelty of our treatment is that the temperature is explicitly put in the 
problem (rather than indirectly through the dependence of the parameters 
of the model). As a consequence we can discuss the critical behavior of the 
transition close to but also far from the critical point. In particular we find a 
softening of the bulk modulus as the transition is approached. 

In many respects the physics is similar to spinodal decompositions of alloys; 
however, an important difference is that in the present case the order param- 
eter is not conserved, slow atomic diffusion is not involved, and the kinetics 
of the problem is very fast. 

The softening effects that we find in a transition that exhibits a discontinu- 
ous volume change, and hence is usually considered as a first-order transition, 
have puzzled the mixed-valence community for decades. They are usually re- 
ferred to as "precursor effects" and often associated with gradually changing 
electronic properties (like "valence"). Sometimes the transition is interpreted 
as a second-order one due to the softening and in contradiction to the obvi- 
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ous discontinuous behavior [31]. Our results provide a simple explanation for 
these precursor effects without invoking peculiar temperature dependence of 
electronic properties and clarify the rather confusing issue of the order of the 
transition. This unusual transition shares characteristics with a second-order 
transition (diverging compressibility) and a first-order transition (discontinu- 
ous order parameter). 

Our analysis has been simplistic, particularly in the fact that we considered 
the two phases to be elastically isotropic and with the same value of the elastic 
coefficients, and the transition to be coherent, namely no structural defects 
(dislocations) were supposed to appear during the transition. Real materials 
are certainly anisotropic and show substantial asymmetry. Although probably 
none of the above assumptions is fully justified, we believe that the features 
obtained are quite robust. Indeed our findings are in good agreement with 
the behavior of real anisotropic mixed- valence systemsfl ,2,27,31]. In these 
systems, the unexpected thermodynamics is a consequence of the interplay 
between strong correlations, determining the double well potential [33, 34], and 
long-range strain effects. It has been suggested that volume instabilities are 
generic for systems with electronic first order transitions [4, 5]. The present ap- 
proach can be seen as a first step to understand this interplay in more complex 
non-isostructural transitions like the magnetoresistant manganites which also 
show volume mismatch among different phases and irreversibility [36]. Inter- 
estingly a volume collapse transition has recently been observed in the parent 
compound [37]. 



7 Appendix 

In order to derive Eq. (4) it is convenient to work with symmetrized strains. For 
d = 3 we define e 2 = (en - e 22 )/V2, e 3 = (en + e 22 - 2633)/^, e 4 = \/2e 23 , 
e 5 = \/2ei3 and e$ = v^2ei 2 . For d = 2 we define e 2 = (en — e 22 )/y / 2 and 
e 3 = y/2e 12 . 

The shear part of the energy can be put as: 

d(d+l)/2 

^ = 7^E E ^kteHO. (io) 

L k i=2 

The St. Venant compatibility condition leads to six equations in d = 3 and to 
one equation in d — 2. For simplicity we restrict to the d — 2 case, the d — 3 
case follows in a similar way. The compatibility condition reads in Fourier 
space: 

C(k) = e(k)-^ - e 2 (k)(k 2 x - k 2 y ) - 2e,^)k x k y = (11) 
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For k = the constrain is automatically satisfied and Hh is minimized by 
e 2 (k = 0) = e 3 (k = 0) = 0. In order to enforce the constraint for k ^ 
we follow Ref. [20,21] and introduce the Lagrange multiplier [20, 21] A(k). 
Minimizing H h + -Mk)C(k) with respect to e 2 and e 3 and eliminating 
A(k) with the aid of Eq. (11) one obtains: e 2 (k) = e(k)(k 2 x - k 2 y )/(y/2k 2 ) 
and e 3 (k) = 2e(k) k x k y / (\/2k 2 ) which allows to eliminate the accommodation 
strains from Eq. (10) to obtain Eq. (4). 
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